Draft version July 19, 2012 

Preprint typeset using I^T^X style cmulatcapj v. 5/2/11 



EXTRACTING PLANET MASS AND ECCENTRICITY FROM TTV DATA 

YORAM LlTHWICK 1 , JlWEI XlE 2 , & YANQIN Wu 2 
Draft version July 19, 2012 

ABSTRACT 

Most planet pairs in the Kepler data that have measured transit time variations (TTV) are near first- 
order mean-motion resonances. We derive analytical formulae for their TTV signals. We separate 
planet eccentricity into free and forced parts, where the forced part is purely due to the planets' 
proximity to resonance. This separation yields simple analytical formulae. The phase of the TTV 
depends sensitively on the presence of free eccentricity: if the free eccentricity vanishes, the TTV will 
be in phase with the longitude of conjunctions. This effect is easily detectable in current TTV data. 
The amplitude of the TTV depends on planet mass and the free eccentricity, and it determines planet 
mass uniquely only when the free eccentricity is sufficiently small. We proceed to analyze the TTV 
signals of six short period Kepler pairs. We find three (Kepler-18,24,25) are consistent with having zero 
TTV phase and are likely devoid of free eccentricities. This result, combined with the observed pile-up 
of Kepler pairs near mean-motion resonances (explainable by resonant repulsion), suggests that the 
orbits of at least some low-mass Kepler planets have experienced substantial dissipation. The fact that 
these pairs likely have zero free eccentricity allows accurate determination of planet masses, subject 
only to uncertainties in transit parameters. The remaining three systems (Kepler-23,28,32) appear to 
have free eccentricities of a few percent. 



1. INTRODUCTION 



The Kepler mission has detected an abu ndance of low- 
mass close-in planets ( Batalha et al.|2 012). Remark ably. 



hundreds of them are me mbers o f plane tary systems ( Lis 
sauer et al.|2011l|Fabrycky et al.|2012bb. These will likely 



prove to be a Rosetta stone for deciphering the dynami 
cal history of planetary systems. 

One of the most intriguing Kepler discoveries is that, 
while the spacing between planets in a system appears to 
be roughly random, there is a distinct pile up of planet 
pairs just wide of certain resonances, and a nearly empty 
gap just narr ow o f them (|Lissauer et al.| 2011| |Fabrycky| 
et al.||2012b[ ). In |Lithwick fc Wu| ()20l2| ), we proposed 
that dissipation is responsi ble for this asymmetry (als o 
see the independent work by Batygin & Morbidelli 2012 ). 
Two planets near resonance repel each other when their 
eccentricities are weakly damped: the inner one moves 
in and the outer one moves out, an effect we termed 
resonant repulsion. The eccentricity of planets near res- 
onances can be separated into two parts: a part that is 
forced by the resonance and is determined by the plan- 
ets' proximity to resonance {forced eccentricity), and a 
part that is unrelated to the resonance {free eccentric- 
ity^ If there is dissipation, it damps away the planets' 
free eccentricities, but the forced eccentricities persist as 
long as the planets remain close to resonance. As the 
dissipation continually acts on these forced eccentrici- 
ties, it extracts energy from the planets' orbits, and in 
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|2000| . But only the resonant contribution plays a role m this pa- 
per. 



doing so pushes apart any planet pair that happens to 
lie ne ar a resonance ( Lithwick fc Wu||2008[ |Papaloizou 
20111. Hence all such planet pairs end up just wide ot 



resonance, naturally explaining the Kepler result. 

If it is indeed resonant repulsion that is responsible for 
the pile up - and to date no other tenable mechanisms 
have been proposed - then there are a number of inter- 
esting implications. First, it implies that before resonant 
repulsion occurred the distribution of spacings was nearly 
uniform, i.e., that planet pairs were placed with little re- 
gard for resonances. Second, it implies that most of the 
Kepler planets suffered a prolonged bout of eccentricity 
damping. Third, it implies that the free eccentricities 
of the planets should be zero today, after the prolonged 
bout of dissipation. This prediction can be tested using 
the transit time variations (TTV) recorded by Kepler, as 
we demonstrate in this paper. 

A transiting planet that has no companion transits at 
perfectly periodic times. But one that has a companion 
deviates slightly from its periodic sche dule because of the 
gra vitational tugs from i ts com panion. Agol et al. ( 2005 \ 
and Holman & Murray (2005) proposed using I T V sig- 
nals to characterize the companions of transiting extra- 
solar planets, and this technique has proven to be highly 
successful both for confirming Kepler candidate s and for 
measuring their masses and eccentricities (e.g. , |Cochran 



et al pilUFord et al.|2012j|Steffen et aL|2012t|FaBrycky 
et al.||2012a ). However, all these studies rely on fitting 
the observed TTV sign als to direct N-body simulations 

(e-g I 
costly 



V eras et al. 2011|) 
Moreover 



Such fits are computationally 
IN -body simulations do not provide a 



dynamically transparent interpretation of the system. 

Here, we focus on near-resonant pairs because the 
nearly coherent interactions in such pairs induce partic- 
ularly large TTV signals. Such pairs account for most 
of the TTV detections to date in the Kepler database. 
Motivated by our earlier work on resonant repulsion, we 
separate the eccentricity into free and forced parts. In- 
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terestingly, in so doing, the expression for the TTV near 
first-order resonance becomes particularly simple. 

This paper is organized as follows. In Section [2] we 
present new analytical formulae for the TTV from two 
near- resonant planets, and show that the results agree 
with N-body simulations. In Section [3j we apply the 
TTV formulae to six planet pairs with published TTV 
data. In Section [4j we discuss our findings and their 
implications. 

2. PLANET PARAMETERS FROM ANALYTICAL TTV 

We consider the TTV signals from two coplanar plan- 
ets that lie near (but not in) a j : j — 1 mean motion 
resonance. The results are derived in the Appendix, and 
summarized in the following. Let St be the deviation 
from exact periodicity in the times of transit of the inner 
planet, and St' that of the outer planet. We have 



V 



(1) 



5t= — e iV + c.c. = \V\ sin(A j + ZV) 
2i 

St' = Xre lX3 + c.c. = \V'\ sin(A J + ZV') , (2) 

where V = \V\e lZV and V = \V'\e lZV ' are the com- 
plex TTV (expression in Eqs. ^n9j), and are nearly con- 
stant; c.c. denotes the complex conjugate of the preced- 
ing term; and 



\i=j\'-(j-l)\, 



(3) 



is the longitude of conjunctions (e.g., Agol et al. 2005), 
where A and A' are the mean longitudes ot the inner and 
outer planet, respectively. At the order of approximation 
to which we work (see Appendix), we may set 

A' = ^(t-T') (4) 



A=|(t-T), 



in Equation ([3]), where P and P' are the periods of the 
inner and outer planet, and T and T' are offsets; all 
four of these parameters are constant. The actual mean 
longitudes differ slightly from the above expressions if the 
planets are eccentric (and that is included in the TTV 
derivation) . 

Equations (JTJ) ([2|) show that TTV signals are sinu- 
soidal, with a period determined by A J . We call this 
the super-period: 



P j 



\j/P> - (j - I)/ P\ 



(5) 



There is a simple geometrical interpretation to A 3 . Since 
A and A' are approximately the angular positions of the 
inner and outer planets, A J gives the angular position 
of both planets whenever they hit conjunction (A = A'). 
Furthermore, X J progresses linearly in time between con- 
junctions. For a near-resonant pair, successive conjunc- 
tions differ only slightly in angular position, and hence 
their super-period is very long. More precisely, we define 
the normalized distance to resonance 



P' 3-1 



P 



1 
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in which case the super-period is 

P' 

P j 



(6) 



(7) 



For example, a planet pair that has a period ratio of 
2.02 is at a distance of A = 0.01 from the 2:1 resonance, 
and its super-period is 50 times longer than the outer 
planet's orbital period. A pair of planets that have A > 
lie wide of resonance, and their A J decreases with time 
(i.e., is retrograde with respect to the orbital motion); 
conversely, a pair with A < lie narrow of resonance 
and their A J is prograde. 

Henceforth it will prove convenient to measure angles 
with respect to the line of sight. In that case, ZV and 
ZV' are the phases of the TTV signals relative to the 
time when the longitude of conjunction points along the 
line of sight (A J =0). In addition, the offsets T and T" 
then have the interpretation of being the time of any par- 
ticular transit of the inner and outer planet, respectively. 

In the Appendix, we derive the expressions for the com- 
plex TTV: 



V = P 



V' = P' 



7rf/3(j_x)i/3A 
\i ( 3 Z, 
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(8) 
(9) 



where /z is the mass ratio of the inner planet to the star 
and fj! that of the outer planet, and / and g are sums 
of Laplace coefficients with order-unity values, as listed 
in Table [3] Note that / < and g > 0. When the pair 
is far from resonance, |A| is order unity, and the TTV 
expression reduces to ~ Pfi. The key new dynamical 
quantity that controls the TTV signal is 



Zs 



free 



/Zfr. 



(10) 



(with Zf ree its complex conjugate), which is a linear com- 
bination of the free complex eccentricities of the two 
planets. We proceed to define and describe the impor- 
tant concept of free eccentricity. 

2.1. Free Eccentricity 
The complex eccentricity of a planet is 

z = ee iro , (11) 

where w is the longitude of periapse. Near a first-order 
mean-motion resonance, this can be decomposed into free 
and forced parts, 



^•frcc T" ^forced • 



(12) 



The planet's forced eccentricity is forced by virtue of its 
companion's proximity to resonance. For the inner and 
outer planets, 



^forced 
''forced 
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(13) 



(Eq. A15|. In the case A > 0, the forced eccentricity of 



the inner planet is aligned with the longitude of conjunc- 
tions, and that of the outer planet is antialigned. The 
situation is reversed for A < 0. 

The free eccentricities can take arbitrary values. They 
represent the degrees of freedom associated with the non- 
circularity of the orbits. Although there are four such 
degrees of freedom (corresponding to e, vj, e',vj'), they 
can only affect the TTV signal through the linear combi- 
nation of Equation (10). The complex free eccentricities 
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Fig. 1. — Schematic of free and forced complex eccentricities of 
a pair of near-resonant planets. The complex eccentricity for each 
planet is a sum of free and forced components. The free component 
is constant (on timescales <C secular times) with arbitrary ampli- 
tude and phase. The forced component has magnitude determined 
by the planets' proximity to resonance, and phase that rotates in 
parallel with the longitude of conjunctions A J , with a super-period 
pi = P'/\jA\. 

precess on the secular timescale, ~ P/ fi, which is much 
longer than the super-period. So during any short time- 
span TTV observation, they can be taken as constant. 
See Fig. [I] for a cartoon illustration of these concepts. 

2.2. Interpreting TTV 

Given two transiting planets, how much can be learnt 
about their parameters from measuring their TTV? The 
analytical expressions allow for straightforward answers, 
sparing one the abstruseness of N-body simulations. If 
both transits are observed, then P, P' , T, T' are known, 
even in the absence of observed TTV. These yield A and 
A J . The TTV signals then allow one to measure four 
quantities, the real and imaginary parts of V and V 
(via Eqs. Urp])- Therefore the four unknown parameters 
(/i, //, ana the real and imaginary parts of Zf ree )could 
in principle be inferred by inverting Equations Q-Q. 
However, as we discuss in the following, degeneracies of- 
ten arise that prevent unique inversion. 

The amplitudes of the complex TTV are approxi- 
mately 



\V\ 

W'\r 



free 



\z. 



free | 



)• 



), 



(14) 
(15) 



after dropping order-unity coefficients. So the TTV am- 
plitude of the inner planet yields an upper limit on the 
mass of the outer planet, // < |A||V|/P, and there is a 
corresponding limit on the inner planet's mass. But the 
value of fjf cannot be extracted from the TTV amplitude 
without knowing the value of |Zf rco |: a smaller // can 
be compensated for by a higher |Zf rC c| without affecting 
the amplitude of V. Similarly, if both planets have mea- 
sured TTV's, one can determine the ratio of their masses 
from the TTV amplitudes (within order unity constants) , 
but not the mass of either individually, without knowing 

\Zf IC c\- 

Can the p/iase.Qbe used to determine mass and free ec- 

4 Recall that we define the phases {ZV and ZV') to be relative 
to the time when the longitude of conjunctions points along the 
line of sight (A J = 0). With this definition, the phases can be 
determined from observed TTV signals with no ambiguity about 
the origin of time or angle. 



centricity uniquely? This is in general impossible. From 
Equations |8])-([9]), the two planets' TTV signals are ex- 
actly out of phase with each other (anti-correlated) in 
either the limit that \Z hcc \ < |A| or |Z frcc | > |A|. This 
feature has been noted before and has been used to help 
confirm some Kepler planets ( [Steffen et al.||2012 ). In ei- 
ther limit, TTV signals only provide three independent 
quantities (|V|, \V'\ and one phase), and hence a degen- 
eracy remains between planet mass and free eccentricity. 
To break it, additional information or assumptions are 
required, e.g., radial velocity measurements or the sta- 
bility of the planet ary system ( |Cochran et al.|[2011l |Fab-| 
2012a|) . An alternative solution would be it 



ktv = Z(V x sgnA) 
b[ tv = Z{V xsgnA) 



rycky et al.||2012a[ ). 

TTV phases can be measured very accurately to discern 
the small deviation from anti- alignment. 

Nonetheless, the phases contain important informa- 
tion, and in certain circumstances they can break the 
degeneracy between mass and eccentricity. For ease of 
discussion, we first define 

(16) 
(17) 

where sgn A = 1 if A > 0, and -1 otherwise. With these 
definitions, cj> ttv = and cj>' ttv — 180° when Zf TCC = 0, 
independent of the sign of A. Note that if 4>ttv = 0, 
then St crosses zero from above to below whenever the 
longitude of conjunctions points along the line of sight 
(regardless of the sign of A); similarly, if <p' ttv = 180°, St' 
crosses zero from below when A J = 0. We consider the 
two possibilities: 

• If the observed TTV's have a phase shift with re- 
spect to A J , that directly implies that free eccen- 
tricities are present. A large phase shift implies 

l^freel > |A|. 

• If there is no phase shift (i.e., 0ttv = arj d 
4>' ttv = 180°), that does not necessarily imply that 
\Zfree\ = 0. Instead, \Z{ TBO \ could be large but the 
phase of Z[ roc vanishes. However, such a coinci- 
dence is unlikely. Even if the phase of Z{ ree van- 
ished initially, secular precession would operate on 
the timescale ~ P/fi to randomize the phase. So if 
the free eccentricity is large (|Zf ree | 3> |A|), TTV 
phases should be randomly distributed between 
and 2n. Conversely, if \Zf Tee \ < |A|, phase shifts 
should be small. If many of the systems observed 
by Kepler have zero or near zero phase shifts, one 
could argue that most of them have small free ec- 
centricities (|-Zfrco]< |A|). Such a result, if found 
true (see Sectional), has interesting implication for 
the dynamical history of these planets. Moreover, 
for the task at hand, it would break the degener- 
acy in the mass determination: when |Zf ree /A| is 



negligible in Eqs. 
determined by the two 



(|9|), fj, and // are directly 
V amplitudes. 



Our above discussion is greatly aided by the analyt- 
ical expressions. Using N-body simulations alone, it is 
difficult to elucidate the degeneracy between mass and 
free eccentricity. Moreover, the fact that TTV phases 
contain important information has hitherto been over- 
looked, only becoming transparent with the analytical 
formulae. 
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We conclude this secti on by comparing with previous 
work. Agol et al. (2005) estimate the TTV amplitudes 
for two near-resonant planets, assuming that the planets' 
initial orbits are circular. In that case, Zf r ee = — Zf OT ced 
initially, and hence Z{ YCC 



(14)-(15) then imply \V\ 



M/A (Eq 
W|A|)(T 



13) 



Equations 



agreement with Equations (29)-(31) of |Agol et al] ( |2005 I 
in the two corresponding limits (\A\ 2 /n <C 1 and ^> 1). 
We note, though, that there is little reason why the ini- 
tial orbits should be circular. If the planets suffer weak 
damping (e.g., by tides or a disk), that would tend to 
damp away the free eccentricities, but the forced eccen- 
tricities would remain intact as long as p roximity to res- 



onance is maintained 
alytically an expression for the T 



Agol et al. (|2005) also derive an- 
V that is valid when 
( see their Append ix). |Nesvorny fc Morbidelli 



(2008) and Nesvorny] ( 2009 ) derived a general analytic 
expression tor the TTV, but their expressions have not 
been applied to Kepler planets due to their algebraic 
complexity, which arises partly because they do not dis- 
tinguish free from forced eccentricities. 

2.3. Testing with N-body Simulations 

We test our TTV expressions with N-body simulations, 
choosing cases that illustrate our discussion above. Fig- 
ure [2] shows transit time variations from two separate 
N-body simulations for a pair of planets that resemble 
Kepler-18 c/d, sh owi ng agreement with the predictions 
of Equations (p~[)— ([To]) . In the top panel the planets have 
zero free eccentricity and hence TTV are in phase with 
\ 3 , while in the bottom panel the inner planet has a 
free eccentricity and hence the signals are not in phase. 
Numerically, we reach the state of zero free eccentricity 
by first weakly damping the planets' velocities to the lo- 
cal circular speed over a few million orbits. Many other 
forms of weak damping will also remove free eccentrici- 
ties. 

Figure [3] illustrates the degeneracies inherent in ex- 
tracting planet parameters from TTV, and how these 
can be partially rem oved with the phase information, as 
discussed in Section [2.2| In the top panel, we perform a 
simulation similar to the one in the top panel of Figure 
[2j with the same periods for the two planets, but in this 
case the masses of the two planets are reduced by factors 
of 0.2 and 0.37. Furthermore, after the free eccentricities 
are damped away the two planets' complex eccentrici- 
ties are increased by Zf rce = z' h . cc = 0.05. The resulting 
TTV's are almost identical to those in the top panel of 
Figure [2] despite the vastly different parameters. Hence 
if one observed such TTV's, one could not determine the 
planets' masses. Nonetheless, the fact that the TTV's 
are in phase with A J in the top panel of Figure [3] is due 
to our judicious choice of Zzf roc = ^ z f ICC = 0. Even 
though they are initially in phase, this cannot remain 
true for long: the black points in the bottom panel show 
the TTV's in the same simulation, but 10 6 days later, by 
which time the two planets' free longitudes of periapse 
have precessed by ~ 70°. As a result, the TTV's are 
no longer in phase with A J . By contrast, the TTV's of 
pairs with zero free eccentricity would always remain in 
phase with the longitude of conjunction (faint green line 
in bottom panel of Fig. 3]). This example shows that if 
a system is observed to have zero TTV phase shift, one 
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Fig. 2. — TTV's of a near-resonant planet pair, without (above) 
and with (below) free eccentricity. Planet periods and masses are 
chosen to be close to those of Kepler-18c/d (see Table 111 with 
A = —0.028 relative to the 2:1 resonance. The black poims show 
the N-body simulated TTV signals for the two planets and the 
green curves the analytic expressions (Eqs. flHToj) . The vertical ar- 
rows show the times at which the longitude ol conjunctions points 
at the observer (X J = 0). The duration between two such arrows 
is the super-period. In the top panel, the pair have zero free ec- 
centricities, so their TTV cross through zero when A J = 0, with 
the inner planet crossing from above, and the outer from below. 
The bottom panel shows the case when the inner planet has free 



eccentricity Zf r , 



0.03i, which is ~ 20 times greater in amplitude 



than its forced eccentricity and is of order |A|. In this case, the 
TTV's are no longer in phase with A J , though the inner and outer 
TTVs are still anti-correlated. The small rapid wiggles in the green 
curves are due to our inclusion of the 3:2 forcing term in addition 
to the 2:1; this clearly has little effect. 



could argue that it is very likely that this system has zero 
free eccentricity (see Fig. [9| , although one cannot be cer- 
tain. But if many systems have vanishing phase shifts, 
it is more probable that the bulk of them have zero free 
eccentricity. One can then determine planet masses from 
TTV data, free from degeneracies. 

3. APPLYING TO KEPLER PAIRS: MASS AND FREE 
ECCENTRICITY 

Equations ([8])~([9]) show clearly how planet mass and 
eccentricity affect the TTV signals. The task to invert 
observed TTV signals to obtain physical parameters is 
now almost trivial. 



terns with confirmed planets (e.g., 


Cochran et al. 


2011 


Fabrycky et al.|2012b Ford et al. 
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2012). Here, we 

-n ' tt 


apply 



ther restrict to systems that have gone through at least 
one complete TTV cycle, leaving us with 6 systems, all 
of which have inner period < 10 days. 

Our results are depicted in Figures |4]-fT0| and Tables 
[T|-[2] All six short period systems appear consistent with 
naving small free eccentricity |^f ree | ^ 0.01, with three 
of them consistent with zero free eccentricity. Before 
discussing these results, we describe our method. 
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TABLE 1 

Complex TTV for Six Kepler Systems 



Kepler ID 




P [d] 


P' [d] 


A 


\V\ [d] 


0ttv 


\V'\ [d] 


0ttv 


18c/d 


2:1 


7.642 


14.86 


-0.028 


0.0037(±7%) 


-4.3° ± 4° 


0.0028(±10%) 


169° ± 5° 


24b/c 


3:2 


8.146 


12.33 


0.0094 


0.028(±20%) 


-3.9° ± 12° 


0.026(±20%) 


180° ± 16° 


25b/c 


2:1 


6.239 


12.72 


0.0195 


0.0026(±20%) 


5.7° ± 12° 


0.00072(±30%) 


200° ± 22° 


23b/c 


3:2 


7.107 


10.74 


0.0077 


0.031 (±40%) 


-68° ± 18° 


0.013(±30%) 


120° ± 17° 


28b/c 


3:2 


5.912 


8.986 


0.013 


0.0082(±20%) 


-50° ± 14° 


0.0086(±30%) 


130° ± 16° 


32b/c 


3:2 


5.901 


8.752 


-0.011 


0.0062(±30%) 


45° ± 15° 


0.0077(±30%) 


228° ± 19° 



Note. — The last four columns are the amplitudes and phases of the complex TT V. e xtra cted from the 
transit times with a least squares fit (see text); the phases are defined in Equations ( |16| -| [T7| l. The phases 
for the first three systems are consistent with zero free eccentricity (flittv = 0, <t>' ttv = lSU ), and those of the 
latter three are not. Errors on the TTV amplitudes and phases are at 68% confidence. 
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Fig. 3. — Degeneracies in extracting planet parameters using 
TTV signals. Here points result from N-body simulations, and 
green curves are theoretical results. The top panel shows a pair 
of planets exhibiting the same TTV signals as those in the top 
panel of Fig. [2] but the ones here have reduced planet masses 
and non-zero free eccentricities (as marked). The TTV's remain 
in phase with A J because we deliberately chose free eccentricities 



with Zzf r , 



■ initially. This panel illustrates that TTV 



measurements alone do not allow unique inference of mass and free 
eccentricity. Additional constraints must be imposed. The points 
in the bottom panel show the same system but 10 6 days later, 
when the phases of the free eccentricities have secularly precessed. 
The TTV's are no longer in phase with A J . The faint green curves 
in this panel are what the TTVs would look like for the system 
with zero free eccentricity (i.e, the system in the top panel of Fig. 
[2| - its TTV remains in phase. Therefore the phase of the TTV 
provides a statistical way to break the degeneracies inherent in the 
TTV's. 

3.1. Extracting the Complex TTV 

Given a sequence of transit times for two planets in 
a system, we obtain the parameters P, P' , T, T' , V, V in 
Equations Q-Q as follows. We write the transit times 
for the inner planet as 

^trans = T + Pitrans + Real(V) sin A- 7 + Imag(y) cos \° 

(18) 

where itrans — 0, 1, ••• is the transit number. Since 
|V| <C P, we proceed in two steps. First, we fit the in- 
ner planet's itrans vs. itrans with a straight line, thereby 
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Fig. 4. — TTV data (black points with errorbars), best-fit theo- 
retical curve, and times when the longitude of conjunction points 
at the observer, i.e., A J = (vertical arrows). These two systems 
have near-zero phase shift, as evidenced by the fact that the inner 
planet's TTV (upper panel) crosses through zero from above at 
times when A J = 0; and the outer planet's TTV crosses through 
zero from below. Hence the se systems likely have zero free eccen- 
tricity. K epler 18 data from ICochran et al.l (12011b , and Kepler 24 
data from |Ford et al.| l |2012[ |. t he former system lies inside the 2:1 
resonance, while the latter one lies outside the 3:2. 
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Fig. 5.— T TV data and fits. S imilar to Figure B] but for Kepler 
25 (data fromjSteffen et al.|2012) ) and Kepler-23 (clada from [Ford| 
|et al,||2012} . Kepler 'Zb nas zero phase shift, but Kepler 23 has a 
non-zero phase shift. 
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Fig . 6.— TTV data and fits fo r KepIer-28 (data from|Steffen et al.| 
|2012| | and Kepler-32 (data from |Fabrycky et al.|2012a| f. Both these 
systems have non-zero phase shirts. 



TABLE 2 
Planet Masses 





analytical Mass [iW®] 


N-body Mass [M ffi ] 


Kepler ID 


^nominal 


^"nominal 


m 


m' 


18c/d 


20.2 ±1.9 


17.4 ± 1.2 


17.3 ± 1.7 


15.8 ± 1.3 


24b/c 


28.4 ±5.9 


50.4 ± 7.9 


56.1 ± 15.8 


102.8 ± 21.4 


25b/c 


7.13 ±2.5 


13.1 ± 2.6 


8.1 ± 3.1 


13.3 ±3.9 


23b/c 


14.7 ±3.8 


55 ± 22 


4.8 ± 15.6 


15.0 ±49.8 


28b/c 


14.8 ±4.2 


22.9 ± 5.6 


3.8 ±6.9 


4.9 ±9.3 


32b/c 


6.0 ± 1.9 


7.59 ± 2.0 


7.2 ±4.1 


5.2 ±3.5 



The nominal masses m nortl i n;1 ] and ml . , 

nominal nominal 



Note. — 

calculated using Equations | |19[ l- |20] l; these likely reflect true 
planet masses for Kepler-18, 24, and 25, but are only upper 
limits for Kepler-23, 28, and 32. For the stellar masses, and 
the N-body deter mined planet mas ses, we use those listed in 
Table 5 of |Fabrycky et al.| |2012a|, except for Kepler-18 for 
which we use values trom Cochran et al. (2011). 



extracting P and T, and also do the same for the outer 



planet, extracting P' , T' 



least squares (e.g. Press et al. 19921 



Our fits are done by linear 
Second^ we fit 



for the four parameters explicit in Equation ( 18 ) (i.e 



T, P, Real(V), Imag(V)) using a second least squares fit 
for the inner planet's transit times; in this fit, X J is de- 
termined by Equations (|3| _ (|4]) at times ttrans, where the 
parameters T,P,T' , P' are the ones obtained from the 
first fits. This refitting for P helps to remove the small 
linear trend that remains after the first fit. We repeat 
for the outer planet. This procedure implicitly assumes 
that the period of the TTV curve is the super-period 
(Eq. [5| . That this assumption is valid for the six Kepler 
pairs of interest can be seen by comparing the fit to the 
data (Figures p}|6]) . The results of the fits are listed in 
TablelU 

To obtain error estimates for the best-fit values, we 
use the covariance matrix provided by the second least- 
square fit and assume that the quoted 1-er errors in the 
transit time data are independent and Gaussian. The er- 
rors in periods (P,P') and fiducial transit times (T, T") 
are negligibly small (fractional error <C 10 -4 ). The 68% 
confidence error ellipses for the values of V and V are 
plotted in Figure [7j and the corresponding 68% error 



0.02 - 



a 

T3 



c 

00 



OA 

CO 




-0.02 - 



-0.02 0.02 

Real(VxsgnA) [day] 

Fig. 7. — Best-fit complex TTVs, plotted in the complex plane. 
Here, sgnA stands for the sign of A - this term is included so 
that pairs with zero free eccentricity would lie on the horizontal 
dotted line with the inner planet to the right, regardless of the 
sign of A. The filled ellipses are the 68% confidence regions for 
the inner planet of each pair (V), and the open ellipses are for the 
outer planet (V), connected here as dumbbells. Kepler 18, 24, 
and 25 are all consistent with having zero free eccentricities, while 
three other systems shown here are discrepant (Kepler 23, 28, 32). 
But even for these, the inner planets still lie mostly to the right 
suggesting that the free eccentricity is small (|-Zf rcc | < |A|). 




-100 100 

0ttv [ de §-] 

Fig. 8. — Phases of the complex TTV. The x-axis is the phase 
of th e inner planet's TTV, and the y-axis is for the outer planet 
(Eqs. |16|]17| . A system with zero free eccentricity would lie at the 
center oi tne plot. The width and height of each rectangle denote 
the 68% confidence limits measured from Kepler data for the 6 
systems in this study. The three curves represent the theoretical 
prediction for the TTV phases when \Z{ ree \ = A (solid curve), 
l-Zfreel = 1.5A (dashed curve) and |Zf rco | = 2A (dot-dashed curve) 
for systems near the 3:2 resonance. It appears that all systems can 
be explained with small free eccentricity, |Zf rco | < |A|. 
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I 1 1 1 I 1 1 1 I 1 1 1 I 
= 0°, UV ' = 180° ±20° 

= 45 °> 0ttv' = 225 ° ±20 ° 
= 70", " = 250° ±20° 




Fig. 9. — The solu tion for planet mass (normalized by the nom- 
inal mass, Eq. \19\ ) versus the magnitude of the weighted free 
eccentricity (normalized by A), for three different measurements 
of TTV phases, for planets near 3:2 resonance. The three ex- 
amples here are chosen to resemble Kepler 24 (blue), Kepler 32 
(orange) and Kepler 23 (light blue). The sizes of the filled circles 
represent the probability of finding the system at this state: at a 
given [Zfreeli we allow the free eccentricity angle to evenly sample 
values between and 27r. The larger circles have more solutions 
that reproduce the observed phases. So at each TTV phase, given 
the uncertainties in phase measurements (taken to be 20° here for 
<t>' v ), there is a range of planet mass and free eccentricity values 
that can reproduce the observations. When </>v ~ 0> the most 



likely solution is Zf TI 



0,m = m nom ; na i, and when 0y 7^ 0, the 



nominal mass is roughly the mass upper limit. 

rectangles for the phases are plotted in Figure [HJ Vanish- 
ing free eccentricity implies that <fr ttv = Z(VsgnA) = 
and (f>' ttv = Z(V'sgnA) — 180°, or equivalently, that St 
crosses zero from above when A J = 0, and St' crosses 
from below. 

We list in Table the masses that would be inferred 
by setting Zf TCC = in Equations (|8])-([9|, i.e., 



^nominal — 



V'A 



nominal " 



P'9 
VA 



Pf 



nf /3 U - 1) 1/3 



(19) 
(20) 



In Figure [9j we argue that these are likely the true planet 
masses for systems showing near-zero TTV phase (Kepler 
18, 24, and 25), whereas for systems with non- vanishing 
phases, they are likely upper limits. Figure [TU] displays 
these nominal masses versus planet radii. 



3.2. 



From Figures [4 
all have very sma 



Systems with negligible TTV phase shift: 
Kepler-18, 24, 25 

-5j we see that Kepler 18, 24, and 25 
phase shifts: St crosses through zero 
from above at the times when A- 7 = 0, and St' crosses 
from below. This can also be seen in Figures [7]-[8j which 
show that at 68% confidence the phases of these systems 
are nearly consistent with zero. We conclude that these 
systems likely have very small free eccentricity, |Zf rcG | <C 
IAI. 



Kepler 18: This system has currently one of the best 
measured TTV, and its TTV phase lies very close to zero. 
If we assume that this system has zero free eccentricity, 
then we may deduce the masses of the two planets with 
small error bars. These values lie close to those obtained 
using an N-body fit (Table [2}. 

But to obtain these mass estimates and small error 
bars, we must assume that the free eccentricities van- 
ish (or \Zfree\ <^ Without this assumption, the 
masses of the planets are degenerate with the free eccen- 
tricity: one could choose an infinite sequence of Zf lee with 
ZZf roc or 7r to reproduce Kepler 18's TTV signals, 
each corresponding to a different set of masses (compare 
the top panels of Figures [2] and [3| . However, the secular 
precession of Zf rcc implies that a specially aligned Zf loc 
would be unlikely. It is far more likely the free eccentric- 
ities are very small, |Zf ree | < |A| ~ 0.03. 

Witho ut explicitly assumin g that the free eccentricity 
is small, Cochran et al. (2011 ) obtain mass estimates and 
errors similar to ours by fitting the TTV signal (alone) 
with N-body simulations (Table [2]). We suggest this re- 
sult is fortuitous: their N-body simulations have not ex- 
haustively searched all possible mass-eccentricity com- 
binations. The same TTV signals can be produced by 
much lower planet masses. In the case of Kepler 18, addi- 
tional information, such as radial velocity measurements 
and theoretical expectations for planet density, allow one 
to eliminate most of the solutions. But this is not pos- 
sible for most systems. This example accentuates the 
importance of an analytical understanding. 

Kepler 25: The TTV phases of Kepler 25b/c are also 
consistent with zero, a strong indication that the planets' 
free eccentricities are very small. From our fits for |V| 
and \V'\ listed in Table [Tl we derive nominal masses for 
the planets as listed in Table [2j These are likely the true 
planet masses. They are consistent with those from N- 
body fits (Table [2]). However, as for the Kepler-18 case, 
the N-body fit without any assumption on the value of 
the free eccentricity is incomplete. An alternative so- 
lution that matches the observations is, e.g., inner and 
outer masses of 0.59M®, 3.0M^ and Zf rcc — —0.05. 

Kepler 24: Similar arguments apply to Kepler 24, 
for which the phases are also consistent with zero. We 
derive masses of 28.4 ± 5.9M e and 50.4 ± 7.9M© for the 
inner and outer plan ets, when we assume th at the free 
eccentricity vanishes. Fabrycky et al. (2012a I list masses 
from N-body fits of 56.1 ± 15.&/W© and 102.8 ± 21.4M e , 
and large eccentricities, of order 0.4 and 0.3 for the inner 
and outer planets. However, the TTV signal suggest the 
absence of free eccentricities 



With planet radii of 2.4 and 2.8i? e ( |Ford et al.||2012 
Batalha et al. 2012[ ), our planet masses imply unrealis- 
tically high densities of ~ 12g/ cm 3 , comparable to the 
density of lead (Fig. 10 1. Although it is possible that we 



5 The solution of |Fabrycky et al.| | |2012afr for Kepler 24 yields 
|Zfree| *C e: their eccentricity vectors nearly cancel in the combi- 
nation fz + gz' . That is likely why their masses are only twice as 
large as ours rather than ~ e/|A| ~ 20 ti mes larger. We note also 
that our TTV super-period di sagrees w ith Fabrycky et al. (2012a i 
by 60% but agrees with that of |Ford et aT]|2TTTIf~fT-onl et al.jfl2uTZ I 
report two more candidates in the system, KU1-1102.04 an d K.U1- 
1102. 03. If real (they are absent from the table in |Batalha et ah] 
|2012| , they are near resonant with the two planets here but with 
Z\ ~ 3% or 4%, respectively. These should not significantly affect 
the TTV amplitude. 
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may be overestimating the masses by assuming Zf rcc = 0, 
we believe a more likely explanation is that the stellar pa- 
rameters we adopt are incorrect. Visual inspection of the 
host light curve of Kepler-24 indicates that it varies at 
the milli-mag level on a timescale of a few days. This 
kind of variability is not expected of solar-type main- 
sequence stars but may be characte ristic of either more 



massive stars (7 Doradus variables, Balona et al. 1994 
Handler|| 1999b , or, perhap s, more evolved stars such as 



giants ( jBalona et al. 2011). The lightcurve also appears 
to be contaminated by nearby sources (G. Handler, pri- 
vate communication). A spectroscopic examination of 
the stellar properties is necessary to ascertain its nature. 
If we assume that the star has a larger radius, i?* = 2i? Q , 
and the same mass, the planets' masses will remain the 
same but their densities will drop by a factor of 8. 

3.3. Systems with TTV phase shift: Kepler-23,28,32 

Three out of six systems we analyze are inconsistent 
with zero free eccentricity: their TTV phases differ sig- 
nificantly from (pttv — and (f>' ttv — 180°. All are near 
3:2 resonances. Planet masses in these cases cannot be 
reliably determined. 

However, as is apparent from Fig. [HJ there is an ab- 
sence of systems with large phase shifts. If these systems 
had l-Zfreel 3> |A|, their TTV phases should be evenly 
distributed between and 360°. The absence of systems 
with large phase shifts suggests that most TTV pairs 
have small free eccentricities, |Zf reG | sa |A| ~ 0.01. This 
absence becomes more striking when we include more 
TTV pairs (Xie et al., in preparation). 

Kepler 23: This pair has the largest phase shift, 
~ 70°. The nominal masses correspond to highplanet 
densities ~ 10g/ cm 3 , comparable to silver. Fig. ^Ishows 
that if one takes |Zf rcc | = 4A = 0.03, then the realmasses 
would be three times lower than the nominal ones. 

Kepler 28 and 32: The nominal masses in these 
two cases all yield planet densities ~ 3g/ cm 3 . The true 
masses are likely smaller (Fig. [9]). However, if the col- 
lective behavior of these TTV pairs indeed implies that 
l-Zfreel ~ |A|, then the true masses lie close to the nomi- 
nal ones. 



4. SUMMARY 



4.1. 



Analytical TTV 

We have derived simple analytical expressions for the 
TTV from two pl anet s near a first order mean motion 
resonance (Eqs. [l} 10). These show that the amplitude 
and phase of the TTV depend on both planet mass and 
free eccentricity. There is an inherent degeneracy be- 
tween mass and free eccentricity which in general pre- 
vents either from being determined independently of the 
other. This degeneracy, however, may be (partially) bro- 
ken under certain circumstances, based on probability 
arguments. 

There is a special moment in time when the longi- 
tude of conjunction points along the line of sight. If 
the phase of TTV is zero relative to this time, then it 
is most likely that the free eccentricity in the system is 
zero. The TTV amplitudes can then be used to uniquely 
determine planet masses. 

When applying this technique to six published systems, 
we find that three of them are consistent with zero TTV 



t 1 — 1 — rn — I 

p[g/crri 3 ]=!=0.1 



0.3 



18d 




J I I I I I l_ 



mass [M e ] 



100 



Fig. 10. — Mass-radius plot for inner and outer planets in the 6 
Kepler systems we examine: The sha ded rectangles show our values 
for the nominal masses (Eqs. |19(]20|, plotted against planet radii, 
encompassing our error estimates, t he nominal masses are likely 
the true masses for Kepler 18, 24, and 25, and upper limits for 
the other systems (as shown by arrows). For comparison, masses 
determined in the literature with N-body simulations are marked 
with filled circles. The high densities for Kepler 24b/c likely stem 
from inco rrect stellar radii (see text). For radii, see references in 
Figs. |4]6| The error estimates on radii are unpublished for Kepler 
24, 25, 23, and 28, and hence assumed to be 10% for this plot. 

phase, while the other three deviate by less than a ra- 
dian. This clustering around zero phase can be most 
naturally explained if all systems have free eccentricity 
of order a percent or less — which is comparable to the 
typical distance to resonance in these systems. For the 
first group, reliable masses can be obtained, while for the 
latter group only mass upper limits. 

Without the analytical TTV expressions and relying 
only on N-body simulations, it would be hard to reach 
these conclusions. 

4.2. Implications of Small Free Eccentricities 

The nearly zero free eccentricities is a qualified suc- 
cess f or the r esonant repulsion theory ( Lithwick & Wu 
20121 see also |Batygin fc Morbidelli|2012p . In that work, 
we explain the observed pile-up near mean-motion reso- 
nances in Kepler planet pairs by invoking damping. We 
found that if Kepler planets have experienced substan- 
tial energy dissipation, but substantially less angular mo- 
mentum damping, the two planets will be repelled from 
each other. This is analogous to the case of an accretion 
disk, where energy dissipation coupled with angular mo 
ment um conservation lead to disk spreading ( Papaloizou 
2011 1. The repulsion effect is particularly prominent for 
planet pairs close to a first-order mean motion resonance, 
leading to a pile-up of planets just wide of resonances. A 
corollary of this theory is that low mass planets in the 
Kepler sample should have little if any free eccentricity. 

Although this appears to be confirmed by three of the 
six systems we analyze, we are puzzled by the small but 
finite free eccentricities in the other three systems. Free 
eccentricity is a record of the dynamical past of plane- 
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tary systems, and the origin of these free eccentricities 
promises further clues about the origin of these planets. 

TTV data have also been reported for Kepler 9b /c and 
Kepler 30b/c ( |Holman et al.|2010|[Fabrycky et al.|2012a| ). 
These pairs contain one or two giant planets. We find 
preliminary evidence that these systems have large TTV 
phases. This, if true, will indicate the presence of large 
free eccentricities in systems of giant planets, in contrast 
to the lower mass planets discussed here. We are cur- 
rently analyzing Kepler public lightcurve data to distill 
more TTV systems. One interesting issue to pursue is 
whether planet pairs at much shorter or much longer or- 
bital periods have the same characteristics as those ana- 
lyzed here. 

Many planet pairs are also near resonance with a third 
planet. This may bring further complications to our 
TTV analysis but is not considered here. 

4.3. Implications of Planet Masses 



All six systems we analyze have orbital periods between 
6 and 15 days, and planet radii r anging from 2 to 7R& . 



We confirm the mass estimates of Cochran et al. (|2011[) 



for K epler 1 8c/d, and the mass estimates of |Fabrycky| 



et al. ( |2012a ) for Kepler 25b/c. Densities of these planets 
range between 0.3 and 2g/ cm 3 . For Kepler 28b/c, 32b/c, 
we obtain mass upper limits that lead to density upper 
limits of ~ 3g/ cm 3 . We also argue that these upper 
limits are likely not too different from the real masses. 

For Kepler 24b/c and 23b/c, the nominal densities are 
~ 10g/cm 3 . For the case of 24b/c, however, we argue 
that the stellar variability likely indicates a larger stellar 
radius. For Kepler 23b/c, further data may reduce the 
error bars in TTV phases and allow a better constraint 
on the value of the free eccentricity. 

Y.L. acknowledges support from NSF grant AST- 
1109776. JWX and YW acknowledge support by NSERC 
and the Ontario government. 
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APPENDIX 
DERIVATION OF ANALYTICAL TTV 

We consider two planets near a first order j : j — 1 mean motion resonance. We assume that the two planets are 
coplanar with each other and with the line of sight. We also assume t hat their eccentricities are small, and we shall 
consider effects up to 0(e 2 ). The total energy, or Hamiltonian, is (e.g., Murray &: Dermo"tt|2000 ) 

GAhm GM^ml Gmm'» , . 

H = — — ri , (Al) 

2a 2a' a' y ' 

where is the stellar mass and the disturbing function due to this resonance is 

R j = fe cos(X j -w) + ge' cos(\ j - w') (A2) 

where the longitude of conjunctions is 

X j ee jX' - (j - 1)A , (A3) 

and we employ superscript j's to refer to the j : j — 1 resonance. In the above, {to, a, e, A, w, to', a', e', A', vo'} are the 
mass and standard orbita l elements for the inner (unprimed) and outer (primed) planets, following the notation of 



Murray & Dcrmott (2000). The coefficients / and g are order-unity, and are functions of j and 

a ee a/a' . (A4) 

These are tabulated for a few resonance in Table |3] 

Given the above Hamiltonian, we proceed to derive and solve the equations of motion in order to obtain the angular 
positions of the two planets (9, 0') as functions of time; these can then be inverted to obtain the times of transit. 

The equations of motion for the inner planet's longitude is 

dX _ 1 2^a dH 

dt ~ y/GM^ m da ^ ' 

^, (A6) 
a 6 

and similarly for A'. We drop the derivative of the disturbing function in the above, and justify this below. Variations 
in the semimajor axes are second order in eccentricity (see below); hence 



X\_( f(t-T)\^( SX 



X' J \ p(t-V) J \5X' 



(A7) 



where the periods are 



The periods will henceforth be treated as constants (i.e., considered to be functions of the unperturbed semimajor 
axes); T and T' are constant reference times; and the variations in the longitudes SX,SX' — 0(e 2 ) remain to be 
determined. 

For the eccentricity equations, we introduce the complex eccentricities 



(A9) 
(A10) 



in terms of which the disturbing function may be written as 



Ri = hfz*+gz'*)e ixi + c.c. , (All) 



2 

where c.c. denotes the complex conjugate of the preceding term. The eccentricity equation for the inner planet is 

dz 1 1% dH 



dt \fGMl m^fa dz* 
to leading order in eccentricity, and similarly for the outer planet, i.e. 



(A12) 



KA-i^fv'f/V^)^ (A13) 



dt \z I P' V fig 
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where the mass ratios are 

/i = m/M*, // = m'/M*. 
Solving the eccentricity equations to first order yields a sum of free and forced terms: 



*free 



1 



A*5 



iA 3 



where the free terms are constant and the normalized distance to resonance is 

. j -IP 1 , 

A EE 1 , 

3 P 

assumed to satisfy |A| <C 1. In the above, we have used the relation 



dt 



2vr 



(jA)-+0(e 2 ) 



(A14) 
(A15) 

(A16) 
(A17) 



We proceed to determine the 0(e 2 ) changes to a, a' , and thereby to obtain 5X, 5X' to this order. To do so, we note that 
the resonant Hamiltonian (Equation ( |Al[ )) has two constants of motion in additio n to the energy: K = A+( j — 1)(T+T') 
and K' = A' — j(T + T'), where the A and T are the usual Poincare momenta ( Murray fc Dermott||2000 ). From the 
constancy of K, the variation in a over the course of the planets' orbits (i.e. 5a) satisfies 



— I" U — I) I m V a ^ + m V a' 



= const , 



(A18) 



discarding terms of higher order in 5a, 8a! , e 2 , e' 2 . Inserting the eccentricities from Equation ( A15 1, only the cross 
terms between the free and forced eccentricities yield time-varying components to 5a, implying 



5a/ a 
5a' j a 1 



v free 'iX 2 



C.C. 



where 



J frcc 



fZf Te e + 5 z free 



is a weighted sum of the two planets' free ecc entr icities. 
The equation for the longitudes (Equation (A6)) becomes, to first order in 8a, 

d ^ 3 2ir 5a 
dt ~ _ 2P"7 ' 
and similarly for the outer planet. The solutions are 



( 5X 
\5X' 



3Z, 



2ijA 



free + c ^ 



Now, to convert from A to 9, we must add the following, valid to first order in eccentricity, 

9 — X = 2e sin(A — zu) 
= -e lX +c.c. 



(A19) 
(A20) 

(A21) 

(A22) 

(A23) 
(A24) 



and similarly for primed quantiti es. S ince we are ultimately interested in transit times, we need only insert for z 
the forced eccentricity (Equation ( |A15 1), because the free eccentricity produces a term with the same period as the 
transits, and hence does produce variations from transit to transit. By the same logic, we may drop the e lX multiplying 
the forced eccentricity if we choose the observer to be at angular position 0. We then have 



- A 

- A' 



1 



+ C.C. 



(A25) 



The final expre ssion for 9,9' to first order is given by the sum of Equation (A7) (after inserting Equation (A22)) with 
Equation (A25), yielding 



P v ; P 1 2i 



2tt 
P 



2tt r V 



P' 2i 



c.c.} 
fee.} 



(A26) 
(A27) 
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TABLE 3 

Coefficients of the disturbing function, defined via Equations {[ATJ-([A2|. Num erica l values expanded to first order in A, 

THE RELATIVE DISTANCE TO RESONANCE (EQ. |A16|) . 





3 ■ .! 


a 


2:1 3:2 4:3 5:4 


-/ 

a 


jb{ /2 + ^Db\ /2 1.190- 2. 20A 2.025 - 6.21A 2.840 - 12.20A 3.650 - 20. 15A 
(j - D&j^ 1 + f + indirect 2: i 0.4284 - 3.69Z^p] 2.484 - 5.99A 3.283 - 11. 9A 4.084 - 19. 86A 


a bi = bs(o) are Laplace coefficients and D = d/dct i 


Murray & Dermott|2000 





|9||-( |10[ l for the 2:1; however, in Equation |8| one should use Zf rec = /zf r 
appropriate for an external perturber. Nonetneless, this distinction is unlikely to be 
interested in the O(A) corrections to the free eccentricities. 

where the amplitudes are 



5cxt 2 f rfln j where g e xt = 0.4284 - 1.17A, 
of practical importance unless one is 



V 



V 1 -- 



P M' -1/2 

7T jA 
P' [i 

IT jA 



I 1 -3/2 3-Zfrcc 

7 2A 



3Z t 



free 



2A 



(A28) 
(A29) 



Note that a = 



((1 + A)^l) 



-2/3 



whence follows Equation (|8|) in the body of the paper. Since we choose 9 = to 

point along tne line ot signt, transits ot tne inner planet occur wnenever v/'Ztt = U, 1, • • • , ana similarly for t he outer 
planet. We conclude that the TTV signals for the inner and outer planets are as given in Equations (|T[)— < 10 ) . 

ASSUMPTIONS AND RANGE OF VALIDITY 

unity. From 



We have assum ed the system is not locked in resonance. Resonance locking occurs when S\ 
Equation (A22|, this implies that our expressions are valid as long as 

efrcc < A 2 //j . 

For typical Kepler systems, |A| > 1% and \i < 10 -4 , and hence this assumption is usually valid. 

In Equation (A6), we dropped the deriv ative of the disturbing function. This term is ~ /xe, and hence is smaller 
than the term that is kept by ~ A (Eq. A21 ). 

We have neglected secular effects. Secular precession occurs on the timescale ~ -P/Mi which is much longer than 
the timescale of the resonant effects considered here, ~ P/A. Hence secular effects only lead to corrections 
to our formulae of order |/i/A| <C 1. Nonetheless, secular precession causes Z{ rcc to precess on thousand-year 
timescales, which randomizes the orientation of Zf rcc . For that reason, we argue in this paper that if the data 
requires ZZf rcc = (or it) for many systems, it suggests that their Z{ roo nearly vanishes. 

The timescale for general relativistic precession is also too long to be important, specifically (duj/dt)^ 1 rj 
Pclc 2 /{6ttGAQ rs 4000 x P 5 5/ d 3 ay year/radian. 

We assume perfect coplanarity. The small inclination dispersion inferred for Kepler multi-planet systems justifies 
this assumption. 



